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Mu-tau neutrino refraction and collective three-flavor transformations in supernovae 
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We study three-flavor collective neutrino transformations in the dense-neutrino region above the 
neutrino sphere of a supernova core. We find that two-flavor conversions driven by the atmospheric 
mass difference and the 13-mixing angle capture the full effect if one neglects the second-order 
difference between the i/ M and v r refractive index. Including this "mu-tau matter term" provides 
a resonance at a density of p « 3 x 10 7 g cm" 3 that typically causes significant modifications of 
the overall v e and v e survival probabilities. This effect is surprisingly sensitive to deviations from 
maximal 23- mixing, being different for each octant. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Neutrinos of different flavor suffer different refraction 
in matter [l|. The energy shift between v e and or 
v T is AV = V2G F Y e n B with Gf the Fermi constant, 
tib the baryon density, and Y e = n e /ns the electron 
fraction. AV is caused by the charged-current i/ e -electron 
interaction that is absent for and v r . For a matter 
density p = 1 g cm -3 we have y/2 Gf^b = 7.6 x 10" 14 eV, 
yet this small energy shift is large enough to be of almost 
universal importance for neutrino oscillation physics. 

In normal matter, p and t leptons appear only as vir- 
tual states in radiative corrections to neutral-current 
and v T scattering, causing a shift AV^ T = V2GfY^ s ub 
between and v T . It has the same effect on neutrino 
dispersion as real r leptons with an abundance Q 
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where n e = n p was assumed. 



2.7 x 1CT 5 , 
(1) 

For the neutron abundance 
= 0.5, but it provides only 



we have used Y n = n n jriB 
a 2.5% correction so that its exact value is irrelevant. 
A large nonstandard contribution to Y° s can arise from 
radiative corrections in supersymmetric models [|[, but 
we will here focus on the standard-model effect alone. 

This "mu-tau matter effect" modifies oscillations if 
AVjjT > Am 2 /2E. For propagation through the Earth 



and for Ami t 



2-3 x 10 eV , this occurs for neutrino 



energies E > 100 TcV. The oscillation length then far ex- 
ceeds rEarth so that AV^ T is irrelevant for the high-energy 
neutrinos that are searched for by neutrino telescopes. 



A "4m/ 2£; implies p w 
10 7 g cm" 3 . Numerical SN density profiles Q re- 



Altcrnatively, the mu-tau matter effect can be im- 
portant at the large densities encountered by neutri- 
nos streaming off a supernova (SN) core [J]. For E = 
20 McV the condition AV^ T 
3 x 

veal that this occurs far beyond the shock-wave radius 
during the accretion phase, but retracts close to the neu- 
trino sphere after the explosion has begun. To illustrate 
this point we show in Fig. [T] the same matter density 
profiles as in Ref. at 1 ms post bounce (red line) 
and at 1 s post bounce (blue line). As a green horizon- 
tal band we indicate the condition AV^ T = Am^ tm /2E 
for a typical range of SN neutrino energies, whereas the 
yellow and light-blue bands indicate the densities cor- 

atm) and 

v r fluxes from a SN are virtually identical, leaving the 
/XT-resonance moot, whereas the H- and L-resonanccs 
cause well-understood consequences that are completely 
described by the energy-dependent swapping probabili- 
ties for v e and v e with some combination v x of the p and 
r flavor Q. Therefore, the traditional view has been that 
genuine three-flavor effects play no role for SN neutrino 
oscillations unless mu and tau neutrinos are produced 
with different fluxes 

In a recent series of papers 0J1, H, EH [Hi El, EE 
E3, El El E3, El, El HE M7\2M it was recognized, 
however, that the traditional picture was not complete: 
neutrino-neutrino interactions cause large collective fla- 
vor transformations in the SN region out to a few 100 km 
(gray shaded region in Fig. [T]). With the exception 
of Refs. [2l|, [23], only two-flavor conversions driven by 



responding to the H-resonance (driven by Am 

l sol 



the L-resonance (driven by Am 2 ol ). The v^, v T , and 



and the small 0i3 have thus far been studied. 
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Am 

We here extend our previous numerical solutions [r7| 
to the case of three neutrino flavors. Our main results 
can be summarized as follows: (i) A two-flavor treat- 
ment indeed captures the full effect if one ignores AV^ r 
and if the ordinary MSW resonances occur outside of the 
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where the Hamiltonian is 2S 
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FIG. 1: Density profiles in terms of the weak potential AV = 
\/~2G-pn e at I ms and 1 s post bounce of the numerical SN 
models described in Ref. [5| (solid lines). The dashed lines 
represent the simplified matter profile of Eq. Q for Ao = 
4 x 10 6 km -1 and A = 5 x 10° km" 1 , used in our numerical 
calculations in Figure [2] As horizontal bands we indicate 
the conditions AV MT = Aml tm /2E, AV = Aml tm /2E, and 
AV" = Am^ ol /2E for a typical range of SN neutrino energies. 
The gray shaded range of radii corresponds to the region of 
collective neutrino transformations. Within the radius r syn 
the collective oscillations are of the synchronized type. 



collective neutrino region, (ii) Including A V^ T strongly 
modifies the v e or D e survival probabilities, influencing 
the neutrino signal from the next galactic SN. (iii) The 
effect depends sensitively on a possible deviation from 
maximal 023- The purpose of our paper is to provide a 
first illustration of these findings that no doubt need to 
be refined in future. 

Our work is organized as follows. In Sec. [ill we present 
the equations of motion which are solved for the three- 
flavor neutrino fluxes in a simplified scenario for the SN 
environment. In Sec. Mil we consider the limit of a van- 
ishing [it matter effect, while our results when it is sig- 
nificant are described in Sec. IIV1 We conclude in Sec. IVl 



II. EQUATIONS OF MOTION 

Mixed neutrinos are described by matrices of density 
p p and p p for each (anti)ncutrino mode. The diagonal 
entries are the usual occupation numbers whereas the off- 
diagonal terms encode phase information. The equations 
of motion (EOMs) are 



Hp = r^p 



-V+V2G1 



d 3 q 

(2tt) 3 



(0q - 0q) (1-Vq-Vp) , (3) 



v p being the velocity. The matrix of vacuum oscillation 
frequencies is f2 p - diag(m 2 , m|, m§)/2|p| in the mass 
basis. The matter effect is represented, in the weak in- 
teraction basis, by V = V2 Gfub diag(F e , 0, Y° s ). For 
antincutrinos the only difference is f2 p — > — tt p . 

In spherical symmetry the EOMs can be expressed as 
a closed set of differential equations along the radial di- 
rection 1171 . |20j . We solve them numerically as previously 



described [17j . now using 3x3 matrices instead of po- 
larization vectors. The factor (1 — v q ■ v p ) in the Hamil- 
tonian implies "multi-angle effects" for neutrinos moving 
on different trajectories [1, H, [ll|. However, for realistic 
SN conditions the modifications are small, allowing for a 
single-angle approximation. We implement this approx- 
imation by launching all neutrinos with 45° relative to 
the radial direction [171 ] . 

As a further simplification we use a monochromatic 
spectrum (E = 20 McV), ignoring the "spectral splits" 
caused by collective oscillation effects [III, E [3 EM H3] • 
Oscillation effects require flavor-dependent flux differ- 
ences. One expects F Vc > F Pe > F Vy . = i 7 ^ = F Vt = F Pt . 
The equal parts of the fluxes drop out of the EOMs, so 
as initial condition we use F v p u ,u r ,p T = and F Uc = 
(1 + e)F Pc with e = 0.25. 

For the neutrino parameters we use Am\ 2 = ^■ ra 1 \ = 



7.6 x 10" 5 eV 



Am?, 



Am 



2.4 x 10" d eV" 



sin O12 = 0.32, sin O13 = 0.01, and a vanishing Dirac 
phase 5 = 0, all consistent with measurements [24|, 
I25I [26|. We consider the entire allowed range 0.35 < 
sin 2 623 < 0.65 because our results depend sensitively on 

e 23 . 

We use a fixed matter profile of the form p oc r~ 3 , 
implying a radial variation of the weak potential of 



AV = Y e X 



R 



(4) 



where R = 10 km is our nominal neutrino-sphere radius 
and Y e = 0.5. In Fig.Q]we show this profile (dashed lines) 
for two different values of Ao = 4 x 10 6 km -1 and Ao = 
5 x 10 9 km -1 . For the former case, the H-resonancc is at 
rn = 1.9 x 10 3 km, the L-rcsonancc at n, = 8.3 x 10 3 km, 



and the /ir-resonance at ry. 



71 km. For the latter 



they are at rn = 2.0 x 10 4 km, = 9.0 x 10 4 km, and 
r^r = 760 km. 1 

The strength of the neutrino-neutrino interaction can 
be parametrized by 



/'o 



= V2G F (F«-F«) 



(5) 



id t Qp = [Hp, g p ] 



(2) 



1 We loosely refer to the radius where Am^ tm /2E = AV^t as the 
\it resonance, although this would be correct only for a small 
vacuum mixing angle in the 23-subsystem. 
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where the fluxes are taken at the neutrino-sphere ra- 
dius R. As in our previous work [l7j we shall assume 
/.to = 7 x 10 5 km -1 . In the single-angle approxima- 
tion where all neutrinos are launched with 45° relative 
to the radial direction [l7| , the radial dependence of the 
neutrino-neutrino interaction strength can be explicitly 
written as 



= Mo 



R 4 



1 



2-R 2 /r 2 



(6) 



While the r~ 4 scaling of /x(r) for r 3> R is generic, the 
overall strength /io depends on the neutrino fluxes and 
on their angular divergence, i.e., on the true radius of the 
neutrino sphere. Our R = 10 km is not meant to rep- 
resent the physical neutrino sphere, it is only a nominal 
radius where we fix the inner boundary condition for our 
calculation. 

The collective neutrino oscillations are of the synchro- 
nized type within the "synchronization radius." For our 
chosen no and for the assumed excess v e flux of 25% we 
find r syn ~ 100 km as indicated in Fig. [T] Collective fla- 
vor transformations occur at r > r syn . Therefore, the fir 
matter effect can be important only if it is sufficiently 
large for r > r syn . 

Figure [1] illustrates that the region where the \it- 
resonance takes place depends on the time after bounce. 
For realistic values of the matter density profile and 
neutrino-neutrino interaction, one expects to lie far 
beyond the collective region at early times. This can be 
inferred from the relative position of r syn and the inter- 
section of the 1 ms profile and the green band. At later 
times though the proto neutron star contracts and r MT 
moves to smaller radii. Eventually r^ T becomes smaller 
than r syn , at which point AV^ T becomes irrelevant. 

In order to mimic these different situations we will us 
a simple power-law matter profile of the form in Eq. ([4]). 
In other words, we will use a mu-tau matter potential of 
the form 



A v = y T eff A 



R 



(7) 



with a fixed F r eff given by Eq. (TTJ) and a variable coef- 
ficient Ao- Therefore early and late times can be repro- 
duced by considering large and small values of Ao, re- 
spectively, as can be seen in Fig. [TJ In other words, we 
will always assume that the ordinary MSW resonances 
are far outside of the collective neutrino region, whereas 
the fir resonance can lie at smaller (vanishing fir matter 
effect) or larger (large /it matter effect) radii than r syn . 



III. VANISHING MU-TAU MATTER EFFECT 

As a first case we consider the traditional assumption 
of a vanishing fir matter effect, which we account for us- 
ing a value of Ao = 4 x 10 6 . We assume an inverted 
Am 2 tm and use a non-maximal value sin 2 623 = 0.4. 



Our numerical calculations for this case are shown in the 
top row of Fig. O The first two panels correspond to 
the radial evolution of the fluxes of the weak interaction 
eigenstates of neutrinos and antineutrinos, respectively, 
whereas in the last two panels we show the evolution of 
the propagation eigenstates. These are the eigenstates 
of flp + V, i.e., of that part of the Hamiltonian Eq. © 
that does not include the neutrino- neutrino interactions. 
In the collective neutrino region, we observe the usual 
pair conversion of the v e and v e fluxes into the /j and r 
flavors. Had we chosen a maximal 23 mixing angle, the 
appearance curves for these flavors would be identical. 

For larger distances the evolution consists of ordinary 
MSW transformations that are best pictured in the basis 
of instantaneous propagation eigenstates in matter (last 
two panels). Beyond the collective transformation region, 
all neutrinos and antineutrinos stay fixed in their propa- 
gation eigenstates. In the weak-interaction basis, on the 
other hand, this implies fast oscillations because we have 
a fixed energy, preventing kincmatical dccohercnce be- 
tween different energy modes. In the panels for neutrino 
and antineutrino interaction states, for radii beyond the 
dense-neutrino region we show as thick lines the average 
evolution as well as the envelopes of the fast-oscillating 
flavor fluxes. 

Another way of describing this evolution is by the level 
crossing schemes of Fig. [3] The upper panel represents 
the case of vanishing AV^ T , corresponding to Fig. 5d of 
Ref. [|| . The central panel represents the case with large 
AVJiT and a 23-mixing angle in the first octant and is 
similar to Fig. 2 of Ref. ^J. In such plots one shows 
the neutrino energy levels as a function of the matter 
density. The continuation of this diagram to negative 
densities gives us the energy levels of antineutrinos: the 
neutrino energy at a negative density really means the 
antineutrino energy at the corresponding positive den- 
sity. For vanishing density (vacuum), we have the three 
vacuum mass eigenstates that are identical for neutrinos 
and antineutrinos. The upper (blue) line corresponds to 
propagation eigenstate 2, the middle (green) line to 1, 
and the bottom (red) line to 3, a scheme representing 
the inverted hierarchy case. These lines represent the 
propagation eigenstates that are adiabatically connected 
for different densities. 

While in vacuum the propagation eigenstates coincide 
with the mass eigenstates, at large densities they cor- 
respond to weak interaction eigenstates. For vanishing 
AVfj, T and at the low energies relevant to our problem, 
the fi and r flavor are not distinguishable so that any 
convenient linear combination can be chosen as interac- 
tion eigenstates. It is convenient to introduce the states 
v 1 and v' r that correspond to a vanishing 23-mixing angle, 
i.e., they diagonalize the 23-subsystem. If the small 13- 
mixing angle were to vanish, the 3-mass eigenstate would 
coincide with v' T . In the upper panel of Fig. [3] and using 



the 



basis, the 2-state connects adiabatically 



to v e and v' , whereas the 3-state connects adiabatically 
to v P and v'. 
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r (km) 

FIG. 2: Radial evolution of the neutrino fluxes, normalized to the initial £> e flux, for a fixed neutrino energy (E v = 20 MeV) 
and an inverted Am^ tm . From left to right: neutrino weak eigenstates, antineutrino weak eigenstates, neutrino propagation 
eigenstates and antineutrino propagation eigenstates. In the first two columns, after bipolar conversions we show the average 
as thick lines and the envelopes of the fast-oscillating curves as thin lines. The top row shows the case of a vanishing /it matter 
effect, while the three bottom rows use a large [it effect with different values for the 23 mixing angle as indicated. 



At the neutrino sphere, the fluxes are prepared in 
i> e and v e eigenstates, which in the case of inverted 
mass hierarchy coincide with the propagation (or mat- 
ter) eigenstates v™ and v™, respectively. In the absence 
of neutrino-neutrino interactions, since the L-resonancc 
is always adiabatic, the v^s leave the star as v%. In the 



case of v e the evolution depends on sin 2 Q 13 @. For 
values larger than 1CT 3 they propagate also adiabati- 
cally (MSW transformation) and escape as D3, whereas 
for values smaller than 10~ 5 the transition at the El- 
resonance is strongly non-adiabatic: there is a jump of 
matter eigenstates from v™ to v™ and the v^s leave the 
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two panels of Fig. [2] that collective pair transformations 
convert the v e and v e fluxes to v' T and v' T as indicated by 
the arrows in the upper panel of Fig. [3] The consequences 
for the subsequent evolution are dramatic. In the case of 
v e a fraction equal to ei*p e stays in and evolves as in 
the absence of neutrino-neutrino interactions, while the 
rest of v e are transformed to ff 1 . As a consequence, the 
final v e flux, normalized to the initial v e one, is expected 
to be approximately e sin 2 812 — 0.08, see thick line in 
the upper left panel in Fig. O In the case of antineutri- 
nos the effect of the collective pair conversion is to inter- 
change the eigenstates in which D e and v' T arrive at the 
H-resonance. Now v e enters the resonance as Df 1 instead 
Therefore, for sin 2 813 > 10~ 3 the resonance is 



of PS 




adiabatic and the v^s leave the star as v\, leading to a 
final normalized flux of approximately cos" 812 — 0.68, 
see the thick line in the second panel in Fig. [2j Instead, 
if sin 2 813 <• 10~ 5 again there is a jump of matter eigen- 
states from to i>™ at the H-resonance. In this case 
v e leaves the star as £>3, leading to a normalized v e flux 
equal to sin 2 613. 

The impact of collective effects is easier to understand 
if we follow the previous literature [ljl EH and observe 
that, in a two-flavor system, the impact of ordinary mat- 
ter can be transformed away by going into a rotating 
reference frame for the polarization vectors. Collective 
conversions proceed in the same way as they would in vac- 
uum, except that the effective mixing angle is reduced. 
Therefore, assuming an inverted hierarchy (IH) for the 
atmospheric mass splitting and a normal hierarchy (NH) 
for the solar splitting, we should consider the level scheme 
as in the upper left panel of Fig. [4] The mass eigenstates 
now approximately coincide with the interaction eigen- 
states because the 23-mixing angle was removed by going 
to the primed states, and the mixing angles involving v e 
are effectively made small by the presence of matter. Of 
course, this level scheme does not adiabatically connect 
to the true vacuum situation. 

The initial state consists of v e and v e and thus es- 
sentially of v\ and v\. Collective conversions driven by 
Am 2 tm then transform v\V\ pairs to ^3^3 pairs in the fa- 
miliar two-flavor way. If both hierarchies are normal, we 
begin in the lowest-lying state and nothing happens. In 
the hypothetical case where both hierarchies are inverted 
(upper right panel in Fig. 0|, we begin in the highest state 
and Am 2 tm drives us directly to the lowest state. Finally, 
if the atmospheric hierarchy is normal and the solar one 
is inverted (lower right panel in Fig. [3J , collective trans- 
formations driven by Am 2 ol take us to the lowest state. 

We have numerically solved the evolution of the three- 
flavor system with a realistic SN matter profile and found 
that the results confirm this simple picture. In a two- 
star as v x . As a consequence, the survival probability is flavor treatmenti the much smaller A mL leads to collec 
P{v e — > ps sin 612 and P{v e — > v e 




FIG. 3: Level crossing scheme of neutrino conversion for the 
inverted hierarchy in a medium with a vanishing AV^ T (up- 
per panel) and a large AV^ T with 23-mixing in the first octant 
(central panel) or the second octant (lower panel). The arrows 
indicate the transitions caused by collective flavor transforma- 
tions. 



sin 813 or 

cos 2 812 for large and small 813, respectively. 

In the presence of neutrino-neutrino interactions, im- 
portant collective effects take place in the inner SN layers, 
where the neutrino density is high. We observe in the first 



sol 

tive transformations at a much larger radius than Am 2 tm . 
In a three- flavor treatment, Am 2 tm therefore acts first 
and takes us directly to the lowest-lying state if the atmo- 
spheric hierarchy is inverted. Otherwise only the hypo- 
thetical case of the lower-right panel in Figgis an exam- 
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FIG. 4: Vacuum level diagram for all hypothetical combina- 
tions of atmospheric and solar mass hierarchies (normal or 
inverted). The 12 and 13 mixing angles are assumed to be 
very small, mimicking the effect of ordinary matter. The ef- 
fect of collective conversions is indicated by an arrow. 



pie where Am^ plays any role. We have numerically ver- 



ified that normal Ar 



combined with inverted Ar 



'sol 

is the only case where A?n 2 ol drives collective transfor- 
mations. Since Am s 2 D | is measured to be normal, the 
previous two- flavor treatments based on Am 2 tm and 0i3 
fortuitously capture the full effect. 

We conclude that in the limit of a vanishing /it effect 
the collective flavor transformations and the subsequent 
MSW evolution factorize and that the collective effects 
are correctly treated in a two-flavor picture. Of course, 
this situation may change if the matter profile is so shal- 
low that the ordinary MSW effects occur in the same 
region as the collective phenomena [2l| . 



IV. LARGE MU-TAU MATTER EFFECT 

Next we calculate the flavor evolution for the same 
model, now including a significant AV Mr , i.e. we assume 
a large Ao- In this case the flavor content of the neutrino 
and antineutrino fluxes emerging from the SN surface 
depend on the strength of AV^ r as well as the choice 
of O23, as can be seen in the corresponding panels of 
Fig. O This dependence is best illustrated with the help 
of the contour plot Fig. [5] where we show the v e and v e 
fluxes emerging from the SN, averaged over fast vacuum 
oscillations. 

If AV^ is so large that the mu-tau effect is strong in 
the region of collective neutrino oscillations, there are two 
stable limiting cases, depending on the 23 mixing angle. 
If the mixing angle is sufficiently non-maximal and in 
the first octant, the collective oscillations transform the 
initially prepared v e and v e fluxes to the propagation 
cigenstates as indicated by the arrows in the middle panel 
of Fig. [31 i.e., we observe pair transformations to v T D T . 
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FIG. 5: Contours in the space of sin 2 623 and Ao for the v e 
(top) and v e (bottom) fluxes emerging from the SN surface 
for both normal (right) and inverted (left) mass hierarchy. All 
fluxes are normalized to the initial v e flux. We show values 
averaged over fast vacuum oscillations. 



This behavior is understood if we assume that in the 
/j,t system we can once more go to a rotating frame and 
now simply imagine that the 23 mixing angle is effectively 
small by the impact of the /jt matter effect. In this case 
1/3 ~ i/ T . Since collective quasi- vacuum oscillations take 
us to the lowest-lying state, the v% state in the inverted 
hierarchy, we are effectively taken to v T v T pairs. Instead, 
if the 23 mixing angle is in the second octant, and v T 
switch roles, explaining that now ^3 « and D3 « v^. 

These are only heuristic explanations. We expect that 
they can be made precise in a true analytic three-flavor 
treatment of collective neutrino oscillations along the 
lines of Ref. [H]. 

For intermediate values of AV^ r and for 23 mixing an- 
gles near maximal, the final fluxes depend sensitively on 
parameters. For intermediate values of AV^ T , there are 
also nontrivial effects for the normal hierarchy. The col- 
lective effects do not place the ensemble into propaga- 
tion eigenstatcs, preventing a simple interpretation. The 
sensitive dependence for intermediate AV^ T is also illus- 
trated in Fig. [6] where we show the emerging average v e 
and v e fluxes as functions of Ao for two values of O23, one 
in the first and the other in the second octant. In Fig. [7] 
we show the same u e and v e fluxes as functions of sin 2 Q23 
for Ao = 1.85 x 10 9 km -1 . One can notice how the fall of 
p ee is not exactly centered at sin 2 #23 = 0.5 but slightly 
shifted to smaller values. This is due to second-order 
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FIG. 6: Fluxes of i/ e (top) and u e (bottom), normalized to 
the initial v e flux, emerging from the SN as a function of Ao 
for a 23 mixing angle in the first (red line) or second (blue 
line) octant. These curves represent cuts through the inverted 
hierarchy contour plots of Fig. [5] at the indicated values of 
sin 2 023- 

corrections to the [lt resonance condition. 

This dependence on the 823 octant leads to a clear 
imprint on the final survival probability. Let us first con- 
sider the first octant. In the case of v e a fraction equal 
to eFp c stays in u^ 1 . However the presence of the \it- 
resonancc in the neutrino channel makes the rest of the 
v e to be transformed to v™. Their subsequent evolution 
would depend on the adiabaticity of the /ir-resonance, 
but it has been shown to be always adiabatic [J- As a 
consequence, the final v e flux is expected to be approxi- 
mately cos 2 0i2 + esin 2 612 ~ 0.76, see thick line in the 
left panel of the second row in Fig. [5] In the case of an- 
tineutrinos the situation is completely analogous to the 
case of vanishing AV^ ir so that P(v e — > v e ) f» cos 2 612 or 
sin 2 813, depending on the value of 813. 

If 823 belongs to the second octant, then the fj,r- 
resonance lies in the antineutrino channel. The crucial 
point is that now all v e are transformed to = be- 
fore reaching the /XT-resonance, see the lower panel in 
Fig. [3] Taking into account that vf does not encounter 
the H-resonance, the survival probability will be always 
P{p e — + v e ) ~ sin 2 812, independently of the value of 
813. 0n the other hand neutrinos do not feel the \it- 
resonancc and therefore their propagation is the same as 
in the vanishing AV^ 1T case. 
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sin 2 9 23 

FIG. 7: Fluxes of (top) and £> e (bottom), normalized to 
the initial Ve. flux, emerging from the SN as a function of 
sin 2 623 for Ao = 1.85 x 10 9 km -1 . These curves represent 
cuts through the inverted hierarchy contour plots of Fig. [5] at 
the indicated value of Ao. 



TABLE I: Summary of the approximate values of the v & sur- 
vival probability for an inverted hierarchy, including or not 
collective effects. Here a small (large) mixing angle O13 stands 
for sin 2 6i3 < 10 -5 (sin 2 613 > 10" 3 ), while a small (large) 
AV Mr represents r MT being smaller (larger) than r syn . 



Collective 
effects 


AV, r 


e 23 


©13 


V, 

leaves as 


P(u e — > U e ) 


no 


any 


any 


small 


VI 


COS 2 ©12 


no 


any 


any 


large 


Vz 


sin 2 6i3 


yes 


small 


any 


small 


i>3 


sin 2 0i3 


yes 


small 


any 


large 


V\ 


COS 2 012 


yes 


large 


< tt/4 


small 


i> 3 


sin 2 0i3 


yes 


large 


< tt/4 


large 


T'\ 


COS 2 012 


yes 


large 


> 7t/4 


any 


i>2 


sin 2 0i2 



We present in TableUa summary of the cases discussed 
so far. One can see the importance of the presence of 
collective neutrino effects, as well as the dependence on 
the strength of the mu-tau matter effect. 

Another interesting feature concerns the position of 
r syn in the presence of a large fir matter effect. As can 
be seen comparing the first two rows of Fig. [21 the ra- 
dius where collective neutrino transformations begin is 
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slightly larger (r s 



115 km) when we include a sig- 



nificant AV^ r . We have checked that, while for a small 
AVj, T the position of r syn is independent of 623, for a 
large AV^ r the onset of bipolar transformations is de- 
layed for nonzero values of ©23. This effect is largest for 
maximal mixing (023 = 7r/4) and symmetric relative to 
023 = 7r/4. A full understanding of this variation pre- 
sumably requires an analytic three-flavor treatment in 
the spirit of Ref . [22j ■ 



V. CONCLUSIONS 

At the relatively low energies relevant for SN neutri- 
nos, charged mu and tau leptons cannot be produced so 
that mu- and tau-flavored neutrinos are not distinguish- 
able in the SN or in detectors. (In the inner core of a 
SN the temperatures may be high enough to produce a 
significant thermal muon density, but this would not af- 
fect the emission from the neutrino sphere.) The impact 
of the small second-order difference between the and 
v T refractive index does not produce observable effects 
as long as one only considers the traditional MSW flavor 
conversion [J]. 

The picture changes if one includes the unavoidable 
effect of collective neutrino transformations in the re- 
gion above the neutrino sphere. If the matter density 
is large enough that AV Mr is comparable to or larger 
than Am^ tm /2£', the survival probability of v e and v e 
can be completely modified and depends sensitively on 
the mixing angle ©23- In future one should also include 
non-monochromatic energy spectra, leading to spectral 
split phenomena that could be more complicated than 
the previously studied two-flavor cases. One should also 
explore the impact of realistic angular distributions and 
of a non-zero Dirac phase in the neutrino mixing matrix. 

Lower-mass progenitors may collapse with a O-Nc-Mg 
core and, on the computer, explode easily because there 
is very little mass in the envelope [27I . Even at core 
bounce and immediately afterward, the density profile is 



so shallow that the ordinary H- and L-resonances may 
occur within the collective neutrino region [2lj . In this 
case the effects discussed here are irrelevant because the 
mu-tau matter effect is negligible. Probably our effects 
are also negligible during the cooling phase of an iron- 
core SN. However, flavor oscillation effects are probably 
largest during the accretion phase of an iron-core SN 
where the flavor dependence of the spectra and fluxes 
is more pronounced than during the cooling phase p8j . 

When it is important, the mu-tau matter effect adds 
one more layer of complication to the already vexed prob- 
lem of collective SN neutrino oscillations. It was previ- 
ously recognized that "ordinary" collective oscillations 
are almost completely insensitive to the smallness of ©13 
as long as it is not exactly zero. Here we have found 
the opposite for the large mixing angle ©23 that is of- 
ten assumed to be maximal. Even small deviations from 
maximal 23-mixing can imprint themselves in the collec- 
tive oscillation effect. Both results are counter-intuitive 
and opposite to ordinary flavor oscillations. 
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